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ABSTRACT 

We present a new determination of the metallicity gradient in M33, based on Keck/LRIS measurements of 
oxygen abundances using the temperature-sensitive emission line [O III] A4363 A in 61 H II regions. These 
data approximately triple the sample of direct oxygen abundances in M33. We find a central abundance of 12 
+ log(0/H) = 8.36 ± 0.04 and a slope of -0.027 ± 0.012 dex kpc" 1 , in agreement with infrared measurements 
of the neon abundance gradient but much shallower than most previous oxygen gradient measurements. There 
is substantial intrinsic scatter of 0.11 dex in the metallicity at any given radius in M33, which imposes a 
fundamental limit on the accuracy of gradient measurements that rely on small samples of objects. We also 
show that the ionization state of neon does not follow the ionization state of oxygen as is commonly assumed, 
suggesting that neon abundance measurements from optical emission lines require careful treatment of the 
ionization corrections. 

Subject headings: galaxies: abundances — galaxies: M33 — galaxies: spiral — HII regions — ISM: abun- 
dances — ISM: evolution 



1. INTRODUCTION 

M33 was one of the first gal axies in whi ch a radial abun- 
dance gradient was recognized (ISearldl971l) . Beginning with 
ISrriitf3?i975l) . a series of studies showed that H II region abun- 
dances in M33 decline by approximately an order of magni- 
tude over the visible extent of the galaxy jKwitter & Allerl 
IT98TI IVflchez et al.l [19881: iGarnett et all fl992h . with a sim- 
ple exponential providing an accurate fit to the gradi- 
ent. A typical value for the gradient fro m this work is 
-0. 12 ± 0.02 de x kpc" 1 (IVflchez et al I fl988l hereafter V88). 
More recently, ICrockett et al.l d2006l hereafter C06) pre- 
sented new metallicity measurements for a small sample 
of H II regions that indicated nearly constant oxygen abun- 
dances across M33, with a best-fit gradient of only -0.012 ± 
0.011 dex kpc" 1 . If confirmed, this very shallow gradient 
would call into question a variety of results concerning, e.g., 
the meta llicity dependenc e of the Cepheid period-luminosity 
relation dLee et al] 12002) and the CO-H2 conversion factor 
(Rosolows ky et al.l l2003) that depend on the abundance gra- 
dient in M33. To complicate matters further, infrared mea- 
surements of the neon abundance gradient, which should be 
less sensitive to systematic effects such as extinction and the 
temperature of the H II region, yield an intermediate gradient 
of -0.034 ±0.015 dex kpc" 1 that is not in very good agree- 
ment with any of the published oxy gen gradients from optical 
data dWillner & Nelson-Patelll2002l. hereafter WNP02). 

In this paper, we present new, direct metallicity measure- 
ments based on detections of the temperature-sensitive [O III] 
A4363 A line in 61 HII regions in M33, approximately 
tripling the sample available in the literature for M33 (or any 
other galaxy). We use these data to provide an improved mea- 

1 Data presented herein were obtained at the W. M. Keck Observatory, 
which is operated as a scientific partnership among the California Institute of 
Technology, the University of California, and the National Aeronautics and 
Space Administration. The Observatory was made possible by the generous 
financial support of the W. M. Keck Foundation. 



surement of the abundance gradient in M33. In §|2] we describe 
the observations, data reduction, and derivation of metallici- 
ties. In O we determine the oxygen abundance gradient in 
M33, attempt to derive neon abundances, and discuss some 
of the implications of our results. We briefly summarize the 
paper in $4] 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Observations 

The M33 Metallicity Project is a long-term program in 
which we plan to obtain deep spectroscopy of ~ 1000 H II re- 
gions in M33 drawn from the unified catalog of Hodg e et"aT] 
(2002). These data will eventually yield a high-resolution, 
two-dimensional metallicity map covering the entire galaxy, 
with a mean sampling rate of > 2 measurements per square 
kpc, and sampling down to < 50 pc in some regions. At 
present, the data set is approximately 20% of its expected final 
size. 

The observations discusse d in this paper were obtained 
with the LRIS spectrograph (lOke et al .ill 9951) on the Keck I 
telescope on 2004 September 22-23 and 2004 October 21. 
Weather conditions during the September nights were clear, 
while the October observations were affected by clouds at 
times. The spectrograph was configured with the 600/4000 
grism on the blue side, a dichroic at 5600 A, and the 900/5500 
and 400/8500 gratings on the red side (in separate exposures). 
The blue-side spectral resolution was ~ 5 A, and spectra typ- 
ically covered 3600 - 5400 A with some variation resulting 
from slit placement on the mask. Slits for which the wave- 
length coverage missed the [O II] A3727 A line have been ex- 
cluded from our analysis. The slits were 1" wide and had 
a minimum length of 10". We used integration times of 15 
minutes per exposure, with 1-4 exposures per slitmask. The 
slits were generally aligned within ~ 10° of the parallactic 
angle during the observations. We observed 1 1 slitmasks in 
the southwest half of the galaxy, with an average of ~ 20 slits 
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on each mask. For the purposes of this paper we use only the 
blue-side data, which include all of the lines ([O II] AA3726, 
3729 A and [O III] AA 4363, 4959, 5007 A) necessary to derive 
oxygen abundances. 

2.2. Data Reduction 

The data were reduced using an IDL pipeline originally de- 
veloped for the Sloan Digital Sky Survey spectroscopic data 
(D. J. Schlegel et al., in preparation) and then modified to deal 
appropriately with multi-slit LRIS spectra. The images were 
bias-subtracted and then flat-fielded using either dome flat ex- 
posures or the internal quartz lamp. We used the flatfield ex- 
posures to identify the location of each slit on the CCD, and 
the remainder of the processing was carried out on each slit 
individually. We produced a two-dimensional wavelength so- 
lution for each slit using exposures of Hg, Cd, Zn, Ne, and Ar 
arc lamps. We then removed cosmic rays from the images and 
coadded the frames together. 

After these basic reduction steps, each slit was examined 
interactively and the targets were identified based on peaks 
in the spatial profile of the H/3 line. We selected an extrac- 
tion region around each target by following the H/3 profile 
down to the background level (which in most cases includes 
emission from the diffuse ionized medium of M33) and set- 
ting the extraction limits where the H/3 intensity reached the 
background value on either side of the peak. A b-spline fit 
to the remainder of the slit produced a 2D sky model cover- 
ing the entire slit. Finally, we subtracted the sky model from 
the data and extracted the H II region spectrum with a boxcar 
algorithm, yielding a ID spectrum and associated uncertain- 
ties. The spectrum was corrected for atmospheric extinction 
and flux-calibrated by comparison to an exposure of the spec - 
trophotometric standard star G191B2B (Ma ssey etaD 1988). 
An example calibrated ID spectrum is shown in Figure Q] 

Our procedure for measuring emission line fluxes in the 
ID spectra began by dividing each spectrum up into a se- 
ries of short (/~ 200 A) spectral windows, each containing at 
least one bright line (except for the regions between 4500 A 
and 4800 A and beyond 5100 A where no such lines exist) 
and performing a constrained linear+Gaussian fit to all of the 
emission lines contained within each window. In each win- 
dow the line separations were fixed at their known values, 
and the widths of the lines were all set equal to the width 
of the brightest line in the window (except for blends such 
as the partially resolved [O II] A3727 A doublet). We then 
measured fluxes by integrating the observed intensity over a 
±1.5 FWHM region centered on each line. Uncertainties in 
the measured line fluxes were calculated by propagating the 
uncertainty associated with each pixel in the one-dimensional 
spectrum. We corrected the measured line fluxes for redden- 
ing and Balmer absorption from th e underlying stellar popu - 
lation using the method outlined bv lOlive & S killman (2001), 
including all of the Balmer emission lines for which we had 
high S/N detections. We used the Galactic reddening law of 
ICardelli et al.l JT989) and assumed that the equivalent widths 
of all of the Balmer absorption lines in a given H II region 
were the same. We experimented with setting the ratios of the 
equivalent widths of the absorption lines equal to their mean 
values averaged over time f or the solar metallicity, continu - 
ous star formation models of Gonzalez Delg ado et al.l (I1999I) . 
but we found that the resulting metallicity changes are always 
less than 0.01 dex. With such minimal changes and with pos- 
sible differences between the stellar populations simulated by 



iGonzalez Delgado et all £l999) and those present in M33 H II 
regions, we proceeded using the simpler model of constant 
Balmer equivalent widths. Corrected line fluxes for each H II 
region with a detection of [O III] A4363 A are given in Ta- 
ble |2] The listed uncertainties in the line fluxes have been 
propagated through the reddening and absorption corrections. 
Balmer line ratios as well as derived reddening and absorption 
values are given in Table [3] 

2.3. Derivation o/H II Region Temperatures and 
Abundances 

We measure oxyge n abundances using the nebular module 
(IShaw & Dufourl 19951) of the STSDAS package in IRAF. We 
model the nebula as consisting of two ionization zones with 
the low- and medium-ionization zones defined by the pres- 
ence of + and O" 1 " 2 , respectively. The temperature of the 
medium ionization zone is calculated using nebular based 
on the standard ratio of the [O III] AA4959,5007 A doublet 
to the [O III] A4363 A auroral line. To determine the tem- 
perature of the l ow-i onization zone, we adopt the scaling of 
ICampbell et all (fl986h : T([0 II]) = 0.70 T([0 III]) + 3000 K. 
We assume that all H II regions are in the low density limit 
(n e < 10 2 cm" 3 ). With these densities and temperatures, we 
measure ionic abundances relative to H + for + and +2 
and assume that all oxygen is found in these two ionization 
states. This analytic approa ch is common to many extragalac- 
tic abundance studies (e.g.. iKennicutt et ai] |2003). but relies 
on several (admittedly well-justified) assumptions. We will 
examine these assumptions in future work with the complete 
sample of H II regions including the full spectral coverage of 
LRIS. 

We restrict our analysis to the 61 H II regions that have de- 
tections of the [O III] A4363 A line with S/N> 3 (more con- 
servative selections do not change our results; see 33. 11 1. The 
locations of these H II regions are superimposed on an Ha 
image of the galaxy in Figure [3] All of our targets from this 
phase of the M33 Metallicity Project are in the southwest por- 
tio n of the galaxy. Pre vious work (with the notable exception 
of Magrini et al.l 12007b has focused primarily on giant H II 
complexes (e.g., NGC 604, NGC 588) in the northeast half 
of the galaxy. The [O III] A4363 A line is detected in - 1/3 
of our targets and this fraction shows no significant variation 
with galactocentric radius. Contaminating supernova rem- 
nants, emission line stars and planetary nebulae are removed 
from the sample based on identification of characteristic spec- 
tral line features and cross-matching with published catalogs 
of these objects. 

As noted in 32.21 uncertainties in the line fluxes are prop- 
agated from single pixel error values. The uncertainties in 
c(H/3) and the equivalent width of the Balmer absorption have 
been included in the errors for the dereddened line fluxes. The 
reported errors in the physical parameters are determined by 
a Monte Carlo propagation of errors through the entire metal- 
licity determination process. The Monte Carlo process is par- 
ticularly well-suited for gauging the widths of skewed uncer- 
tainty distributions in the derived properties from low signal- 
to-noise lines. For each H II region, we generate 100 sets 
of trial line fluxes by adding a normal deviate times the flux 
uncertainty to each line flux and re-analyzing the spectrum 
using nebular. We also include a random fluctuation in the 
temperatur e of the low ion i zation zone derived from the rela- 
tionship of |Ca^be|retal]([T98a). Ba sed on the magn i tude o f 
the uncertainties found in the work of lKennicutt et alj ([2003 ), 
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FIG. 1 . — Example spectrum of H II region 60 in Table|2](BCLMP 266). This spectrum has the median S/N (8.4) for the [O III] A4363 A line in our sample. 
F\ is given in units of 10~ 16 erg s -1 cm -2 A -1 . 



the lcr error in T([0 II]) is 300 K. The final uncertainties on 
the abundances are defined as the lcr widths of the Monte 
Carlo abundance distributions (see also C06). For reasonably 
strong [O III] A4363 A emission (S/N> 10), precisions in the 
derived metallicity are < 0.08 dex. For the four H II regions 
in our sample that have been studied by previous observers, 
the abundance determinations agree to within the uncertain- 
ties (see TableQ]). 

2.4. He II Zones 

We detect He II A4686 A from 5 of the 61 H II regions in 
our sample at a typical level of 2.5% of the H/3 flux (objects 
3, 8, 16, 18, and 29 in Tabled- The emission is almost cer- 
tainly nebular as it is spectrally unresolved but spatially re- 
solved in these sources (C06). The presence of He II emission 
indicates that the simple two ionization zone model adopted 
in the metallicity estimates is somewhat inaccurate. A high 
ionization region implies the presence of O" 1 " 3 , which is not 
accounted for in our abundance determination. Consequently, 
the derived oxygen abundances will be systematically low by 
an unknown amount. To estimate the magnitude of this ef- 
fect, we calculate He +2 /He + based on the ratio of the He II 
A4686 A and He I A4471 A lines. This ratio is only approxi- 
mate because the He I line is likely contaminated by stellar ab- 
sorption. The derived ratio is roughly equal to 3+ /(0 + +0 2+ ), 
providing an estimate of the contribution of 3+ . We therefore 
expect that we have underestimated the oxygen abundance by 
~ 0.01-0.06 dex in these regions. Since we are not sure of 
the influence of stellar absorption lines, we add (in quadra- 
ture) this derived increase to the reported oxygen uncertainty 
for these five regions. 

3. RESULTS 

3.1. The Oxygen Gradient 

Figure |2] shows the radial distribution of gas phase oxygen 
abundances in M33. We have assumed a t hin disk geome- 
try ba sed on a distance to M3 3 of 840 kpc ( Freedma n et al.l 
l200l . i = 52° and PA = 22° (ICorbelh & Saluccill2000l) .~An 
ordinary, linear least-squares fit to the data, weighting each 
point according to its inverse variance, produces a gradient of 
12 + log(0/H) = (8.38 ±0.03) -(0.032 ± 0.006) R kpc , but the 
residuals from this fit are much larger than would be expected 

from the measurement uncertainties (\ 2 = 3.2). Since it is ev- 
ident that an exponential decline with radius is a reasonable 
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FIG. 2. — Abundances of 61 Ft II regions in M33 as a function of galac- 
tocentric radius. A linear gradient with a slope of —0.027 dex kpc~' is fit to 
the data (solid line). Regions with significant He II A4686A emission are 
indicated with open symbols. 

functional form to describe the data, the high x 2 value indi- 
cates the presence of substantial intrinsic scatter around the 
relation. 

In this situation, accurately es timating the fit parameters 
must be done using the method of lAkritas & Bershadyl (1 19961 
hereafter AB96). We perform a linear regression with the 
AB96 technique for weighted least-squares fitting in the pres- 
ence of intrinsic scatter in the data. The AB96 technique es- 
timates the intrinsic variance and adds this variance to that of 
the data to produce the appropriate weights for the fit. Since 
the intrinsic variance is larger than the statistical uncertain- 
ties for the M33 abundance data, there is less variation in the 
weights than in the measurement uncertainties. Even though 
[O III] A4363 A detections with S/N as low as 3 result in rela- 
tively large uncertainties in the derived abundances, including 
these data improves the quality of the fit by increasing the 
sample size. This technique yields an abundance gradient of 

12 + log(O/H) = (8.36±0.04)-(0.027±0.012)fl i7X .. (1) 

The implied intrinsic variance in the population (i.e., the 
scatter in the metallicities at any given radius) is 0.11 dex, 
larger than the precision of most of the measurements. We 
check the uncertainties on the central abundance and the 
slope of the abu ndance gradient by bootstrapping the data 
dPress et al.lll992l) . which yields uncertainties < 10% higher 
than those derived from the AB96 formalism. Although all 
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TABLE 1 

Comparison with Literature Metallicity Measurements 



H II region KA81 V88 C06 M07 Our 

Abundance Abundance Abundance Abundance Abundance 



NGC 588 ••• 8.30 ±0.06 ■•• ••• 8.315 ±0.061 

MA 2 8.38 8.44 ±0.15 ■■• ••• 8.334 ±0.083 

BCLMP90 ••• ■•■ 8.50 ±0.06 ••• 8.503 ± 0.057 3 

BCLMP 218 — 8.25 ±0.05 8.157 ±0.059 

REFERENCES. — KA81 = Kwitter & AHei 1981; V88 = Vflchez et al. 1988; C06 = 
Crockett et al. 2006]; M07 = IMagrini et aU2007l 

a Weighted average of two independent measurements at the same position; see Table|2] 



the included data appear to have a statistically (and visibly) 
significant [O III] A4363 A detection, performing the analysis 
only on the 23 H II regions with S/N> 10 in this line yields an 
indistinguishable gradient of 12 + log(0/H) = (8.41 ±0.07)- 
(0.039 ±0.021) R kpc . We note that all of our targets are lo- 
cated in southwest half of the galaxy (see Figure[3]l, and there- 
fore that the other half of the galaxy could conceivably have 
a different gradient, but we regard this possibility as unlikely. 
Finally, we have correlated our derived metallicities and their 
residuals around a radial gradient against numerous param- 
et ers in our analysis such as L(Ha) taken from the catalog 
of lHodge etalJ (120021) . c(H/3), the equivalent widths of the 
Balmer absorption and H(3 emission, and the oxygen ioniza- 
tion correction factor. We find no correlation with the residu- 
als that would indicate systematic errors in our methods. 

3.2. Complications in Measuring the Neon Gradient 

The slope of the neon abundance gradient is expected to be 
equal to that of oxygen since these elements are produced in 
nearly equal proportions for all channels of stellar nucleosyn- 
thesis. WNP02 derived neon abundances from Infrared Space 
Observatory spectra and found a gradient consistent with our 
oxygen data. Optical observers commonly derive a neon gra- 
dient based on the [Ne III] A3869 A line, yielding an Ne +2 
abundance. The fraction of neon found in the Ne +2 state is of- 
ten assumed to be the same as the fraction of oxygen in +2 , 
so that an identical ioniza tion correction factor (ICF) can be 
used for the two elements (IStasinska et al.ll200ll C06). Using 
the simple ICF results, we found a significant gradient in the 
Ne/O ratio, contrary to expectations. We find cause for con- 
cern with this assumption for the neon ICF when we derive 
the abundance of Ne + - using the standard method and plot the 
Ne +2 to O" 1 " 2 ratio as a function of the fraction of oxygen found 
in the doubly ionized state (FigureE}. This ratio should be in- 
variant with oxygen ionization fraction and equal to the Ne/O 
ratio, yet the plot shows a significant slope, with Ne +2 /0 +2 
varying by a factor of ~ 4 over the observed range. Since the 
relationship between the two variables is not known a priori, 
we fit a linear regression using the BCES method of AB96 to 
account for the uncertainties in both directions. The resulting 
fit gives: 

log(Ne +2 /0 +2 ) = (-1 .45 ± 0.10) + (1.27 ± 0.22) +2 /0. (2) 

This analysis finds that the slope of the relationship is strongly 
inconsistent with the value of zero expected from using the 
same ICF for neon and oxygen. The fraction of doubly- 
ionized oxygen is a proxy for the ionization state of the neb- 
ula, so the trend indicates less Ne +2 relative to O" 1 " 2 in low ion- 
ization nebulae. This is the expected direction for this trend 
since the third ionization potential of neon (41 eV) is larger 
than that of oxygen (35 e V). More recent work bv llzotov et al.l 



d2004t 120061) suggests that applying the oxygen ICF to neon 
may not be valid in lower ionization regions because of charge 
transfer reactions between 2+ and atomic hydrogen. These 
discrepancies could also be caused by a failure to account ap- 
propriately for photons with hv ~ 40 eV in the O-star atmo- 
s pheres used in photoionizatio n codes (C06). 

iPerez-Montero et ail (|2007) propose a nonlinear ICF based 
on their own photoionization models: 

ICF(Ne 2+ ) = 0.753 + 0.142jc+^i, (3) 

x 

where x = 2+ /(0 + + 2+ ). Adopting this ICF, we recalcu- 
late neon abundances and find no significant gradient in Ne/O 
with radius. For our sample of H II regions, we measure 
(log(Ne/Q)) = -0.69 ±0.0 7, consistent with the results of 
IPerez-Montero et ail (120071) for giant extragalactic H II re- 
gions. We conclude that using the same ICF for neon and 
oxygen is inappropriate for the H II regions we observed and 
that more sophisticated ICFs seem to produce reasonable re- 
sults for generic H II regions in M33. 

3.3. Comparison to Previous Gradient Measurements 

Our measured gradient of -0.027 dex kpc" 1 is consistent 
within the uncertainties with the neon gradient of WNP02 
(-0.034 ± 0.015 dex kpc" 1 ) and marginally so with the oxy- 
gen gradient of C06 (-0.012 ± 0.01 1 dex kpc" 1 ). However, it 
is not consistent with that of V88 (-0.05 ± 0.01 dex kpc" 1 for 
the outer regions and -0.10 dex kpc" 1 overall after rescaling 
to our assumed distance) and various older studies. We note 
that the steep inner gradient measured by V88 is driven by 
abundances determined from photoionization models, which 
may be systematically high compared to electron temperature 
abundances. The apparent discrepancy in abundance gradient 
likely results, at least in part, from an underestimation of the 
uncertainties imp arted by intrinsic scatter in the H II region 
abundances (see 83.41 i. The intrinsic scatter in our optical H II 
region abundances (0. 1 1 dex) is comparable to the value of 
0.07 -0.10 dex reported by WNP02. 

WNP02 suggested that a linear gradient may not be an 
appropriate description for the enrichment of M33, arguing 
for a flat distribution of abundances with a step down at a 
galactocentric radius of 4 kpc. We fit a linear gradient to 
the 38 H II regions inside R ga \ = 4 kpc and find a slope of 
-0. 015 ± 0.024 dex kpc" 1 , which is consistent with no gradi- 
ent. IStasinska et al.l £2006) also claimed an abundance plateau 
in the inner galaxy based on several sub-solar metallicity plan- 
etary nebulae in the region. We measure a mean metallicity in 
the inner galaxy (R gSL i < 4 kpc) of 8.30; in the outer galaxy, the 
mean is 8.21, implying a jump of the same order as, but some- 
what smaller than, that of WNP02 if a step function metallic- 
ity profile holds. 
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FIG. 3. — Locations of the 61 H II regions with oxygen abundance determina tions superimposed on a continuum-subtracted Ha image of the southwest portion 
of M33. The image is produced from the survey data of Massey et al. 1 2006). Contours of constant galactocentric radius are shown as dotted ellipses with a 
spacing of 1 kpc. 
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FIG. 4. — Ratio of doubly-ionized neon to doubly-ionized oxygen plotted 
as a function of +2 ionization fraction. Based on photoionization models, 
this value should be constant and equal to the Ne/O ratio. However, there is 
a clear trend of the ratio between the doubly ionized species and the degree 
of ionization in real H II regions. The dashed line shows the fit given in 
Equation ff] Regions with significant He II A4686A emission are indicated 
with open symbols. 



There now appears to be a growing consensus that the oxy- 
gen abundance gradient in H II regions in M33 is relatively 



shallow, in conflict with the decades-long assumption of a 
gradient at least as steep as -0.10 dex kpc -1 . If we ignore 
the innermost H II region observed by V88, which does not 
have a detection of [O III] A4363 A or any other temperature- 
sensitive lines, then all of the modern H II region studies re- 
port gradients of - 0.06 dex kpc" 1 or shallower (V88, C06, 
iMagrini et alj2007l) . We show in ^3~4l that the results of these 
studies, which include 6- 14 H II region abundances, are con- 
sistent with our measurements once their small sample sizes 
are taken into account. 

Other types of sources besides H II regions are increasingly 
being employed to study abundances in M33. lUrbaneia et alj 
(2005) used quantitative spectroscopy of OB supergiants to 
determine an oxygen gradient of -0.06 ± 0.02 dex kpc" 1 . 
Since these stars are very young, they are expected to trace 
the interstellar medium (ISM) abundances. If the OB stars 
have a similar intrinsic abundance scatter to the H II regions 
we measu re, then this result w ould be compatible with our 
findings. Magr inTet alJ ({2007) argue for a steepening of the 
gradient at small radii based o n the young star met allicities 
in the inner 2 kpc of the galaxy (lUrbaneia et al.ll2005l) . but we 
find no evidence for such a stee pening based on the H II region 
data alone. Planetary nebulae (Ma grini et al.l 2004) still yield 
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FIG. 5. — Distribution of the abundance gradients that would be mea- 
sured for samples of 10 H II regions drawn randomly from our sample of 
61. The large width of distribution shows the uncertainty imprinted upon the 
abundance de t ermin ations by the underlying v arianc e. The measurements of 
IVilchez et ai] 0.988; V88) and lCrocketteTai] <2006l C06) are indicated and 
are consistent with the distribution given their small sample sizes. 

a steep oxygen gradient (-0.14 dex kpc" 1 ), but this result is 
dependent on a single nebula at large radius; the remaining 
objects clearly have a much shallower gradient. Finally, su- 
pernova remnants have also shown a shallow oxygen gradient, 
although the gradients i n other elements appear to be steeper 
dBlair & Kirshnerjfl98l . 

3.4. Intrinsic Scatter and the Determination of Gradients 

Even after accounting for uncertainties in the stellar ab- 
sorption and reddening corrections, we find that there is an 
intrinsic scatter of 0.11 dex around the gradient that is un- 
explained by the measurement uncertainties. Such scatter is 
comm only seen in gradient d eterminations (e.g., in the Milky 
Way; lAfflerbac h et al.l 1 1997b and may result from metallic- 
ity fluctuations in the interstellar medium. Regardless of its 
source, gradient determinations made in the face of signifi- 
cant scatter coupled with a limited number of observations 
may produce widely varying results. The historical evolution 
of the gradient determination in M33 (one of the best-studied 
galaxies), ranging over nearly an order of magnitude, serves 
as a cautionary example. 

Only large numbers of measurements can overcome the un- 
certainties engendered by the intrinsic variance and relatively 
shallow gradient in M33. In Figure [5] we show the likeli- 
hood of measuring a particular value of the gradient if only 
10 H II regions from our sample were used in the determina- 
tion of the gradient. We model the behavior for 10 objects 
since this is a typical number of targets studied in estimates of 



galactic abundance gradients, such as those of V88 and C06. 
The broad distribution of the derived gradients (standard de- 
viation of 0.03 dex kpc" 1 ) suggests that uncertainties in the 
gradients are systematically underreported in galaxies with 
significant scatter around their gradients. Even with a sam- 
ple several times larger than previous studies, we can only re- 
port the presence of a gradient in M33 with 23a confidence! 
Other galaxies may have better established gradients owing to 
smaller intri nsic scatter and/or a s teeper abundance gradient 
(e.g.. MIOh lKennicutt et alJl200l . 

4. SUMMARY AND CONCLUSIONS 

Using data from the M33 Metallicity Project, we have ob- 
tained direct oxygen abundance measurements based on de- 
tections of the [O III] A4363 A line for 61 H II regions in the 
southwest half of M33. This data set approximately triples the 
sample of M33 H II region metallicities in the literature and 
makes M33 by far the best-studied external galaxy in terms of 
its ISM abundances. We presented a refined oxygen gradient 
for M33 and found an exponential abundance profile with a 
gradient of -0.027 ±0.012 dex kpc" 1 and a central abundance 
of 12 + log(0/H) = 8.36 ±0.04. There is an intrinsic scatter 
of 0. 1 1 dex around the trend, which complicates the measure- 
ment of the shallow gradient in the galaxy. When this scatter 
is accounted for, previous studies (which generally derived 
steeper gradients) were shown to be consistent with our mea- 
surements. We also noted that there is a significant correlation 
of the Ne +2 /0 +2 ratio with the fraction of oxygen found in the 
doubly ionized state, which suggests that neon determinations 
based solely on optical lines require more sophisticated ion- 
ization correction factors than are typically assumed. 
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TABLE 2 
H II Region Line Fluxes 



Obj. 


Coordinates" 


Name b 




[O II] C 


[Ne III] C 


[O III] C 


[O III] C 


[O III] C 


T e ([0 III]) 


12+log([0/H]) 




(^2000:^2000) 




(kpc) 


A3727 A 


A3869 A 


A4363 A 


A4959 A 


A5007 A 


(K) 




1 


01:33:48.4 +30:39:36 


B0043b 


0.19 


1.42±0.01 


0.086 ± 0.002 


0.010 ±0.001 


0.680 ± 0.002 


2.039 ± 0.006 


9200 ± 300 


8.214 ±0.087 


2 


01:33:47.8 +30:38:38 


B0029 


0.30 


1.58±0.01 


0.0437 ± 0.0007 


0.0046 ± 0.0005 


0.409 ±0.001 


1.229 ±0.003 


8600 ± 300 


8.211 ±0.097 


3 


01:33:44.0 +30:38:49 


B0038b 


0.53 


0.85 ±0.03 


0.51 ±0.02 


0.045 ±0.004 


2.25 ±0.03 


6.50 ±0.07 


10200 ±400 


8.391 ±0.090 


4 


01:33:50.1 +30:37:30 


BOO 16 


0.56 


2.29 ±0.02 


0.023 ± 0.001 


0.005 ±0.001 


0.231 ±0.001 


0.694 ± 0.003 


10300 ±900 


7.98 ±0.19 


5 


01:34:00.3 +30:40:48 


B0079c 


0.73 


1.26 ±0.05 


0.098 ± 0.004 


0.012 ±0.001 


1.02±0.01 


3.05 ±0.03 


8700 ± 200 


8.386 ±0.092 


6 


ai to a c c , in t/" Art 

01:33:45.5 +30:36:49 


B0027D 


0.74 


2.75 ±0.01 


n mm 1 n r\f\f\s 

0.0303 ± 0.0006 


0.0032 ± 0.0004 


f\ on A I f\ f\r\f\c> 

0.3214 ±0.0008 


n nc"7 1 n nno 

0.957 ± 0.002 


o a nn i i nn 

8400 ± 300 


o o"7 i n 1 1 

8.37 ±0.11 


7 


01:33:34.9 +30:37:05 


T~i nn n o 1 

B0033b 


1.25 


l.80l ±0.009 


0.0678 ± 0.0007 


r» f\f\C A 1 f\ f\i\f\ A 

0.0054 ± 0.0004 


n con 1 f\ nn 1 

0.589 ±0.001 


i ti 1 i n nm 

1.771 ±0.003 


o o nn 1 o nn 

8200 ± 200 


o t nn i n non 

8.399 ±0.080 


8 d 


01:34:04.2 +30:38:09 


B0090 


1.28 


l.25±0.05 


0.57 ±0.02 


0.049 ± 0.003 


2.55 ±0.02 


7.73 ± 0.07 


10000 ±200 


8.506 ±0.070 


8 d 


01:34:04.2+30:38:09 


B0090 


1.28 


l.8l ±0.03 


0.652 ±0.009 


0.048 ±0.002 


2.47 ±0.02 


7.39 ±0.06 


10100 ±100 


8.500 ±0.064 


9 


01:34:04.3 +30:38:07 


B0090 


1.30 


2.7 ±0.2 


0.69 ±0.05 


0.041 ±0.009 


2.56 ±0.07 


7.4 ±0.2 


9600 ±500 


8.62 ±0.14 


10 


01:34:07.9 +30:39:14 


C204b 


1.43 


2.7±0.l 


0.072 ±0.004 


0.010 ±0.003 


0.602 ±0.007 


1.83 ±0.02 


9500 ±600 


8.27 ±0.15 


11 


01:33:47.7 +30:33:38 


BOO 18 


1.56 


l.35±0.04 


0.220 ±0.007 


0.014 ±0.004 


0.88 ±0.01 


2.54 ±0.03 


9600 ±800 


8.19±0.17 


12 d 


01:33:42.9 +30:33:30 


B0023 


1.57 


l. 63 ±0.03 


0.477 ±0.009 


0.021 ±0.002 


1.725 ±0.009 


4.93 ±0.02 


8800 ±300 


8.55 ±0.11 


12 d 


01:33:42.9 +30:33:30 


B0023 


1.57 


l. 81 ±0.04 


0.47 ±0.01 


0.031 ±0.004 


1.80 ±0.01 


5.30 ±0.03 


9700 ±300 


8.444 ±0.081 


13 


01:34:01.8 +30:35:49 


B0001 


1.57 


1.68 ±0.02 


0.035 ±0.002 


0.008 ±0.002 


0.379 ±0.003 


1.071 ±0.005 


10400 ±700 


7.93 ±0.14 


14 


01:33:50.3 +30:33:42 


BOO 15b 


1.59 


3.2±0.2 


0.68 ±0.05 


0.023 ±0.007 


2.01 ±0.04 


5.9±0.1 


8600 ±800 


8.73 ±0.23 


15 


01:34:10.5 +30:39:16 


C208 


1.64 


2.89 ±0.02 


0.0167 ±0.0008 


0.0063 ±0.0006 


0.294 ±0.001 


0.947 ±0.002 


10200 ±300 


8.101 ±0.083 


16 


01:34:00.3 +30:34:17 


COOlAb 


1.85 


1.22 ±0.02 


0.63 ±0.01 


0.041 ±0.003 


2.64 ±0.02 


8.03 ±0.06 


9400 ±200 


8.610 ±0.080 


17 


01:33:54.1 +30:33:10 


BOO 13c 


1.86 


1.68 ±0.01 


0.0394 ± 0.0009 


0.0049 ± 0.0008 


0.451 ±0.002 


1.349 ±0.004 


8500 ± 400 


8.26±0.12 


18 


01:33:34.1 +30:32:13 


B0208f 


2.02 


2.69 ±0.02 


0.127 ±0.002 


0.0147 ±0.0009 


0.580 ±0.002 


1.742 ±0.005 


10900 ±200 


8.071 ±0.063 


19 


01:34:16.8 +30:39:21 


B0723 


2.15 


2.4 ±0.2 


0.088 ±0.009 


0.014 ±0.003 


0.78 ±0.02 


2.30 ±0.05 


9800 ± 700 


8.24±0.16 


20 


01:33:27.7 +30:31:56 


B0213b 


2.30 


4.10±0.06 


0.493 ± 0.008 


0.034 ±0.002 


1.718 ±0.009 


5.13 ±0.03 


10100 ±200 


8.495 ±0.062 


21 


01:34:15.6 +30:37:11 


MA2 


2.32 


1.332 ±0.005 


0.0626 ± 0.0005 


0.0053 ± 0.0003 


0.5695 ± 0.0009 


1.733 ±0.002 


8200 ± 100 


8.334 ±0.083 


22 


01:34:09.5 +30:34:18 


C184A 


2.42 


2.31 ±0.04 


0.049 ± 0.003 


0.008 ± 0.002 


0.470 ± 0.004 


1.401 ±0.009 


9600 ± 800 


8.17±0.18 


23 


01:33:23.4 +30:31:35 


CI 29a 


2.55 


2.65 ±0.01 


0.0859 ± 0.0009 


0.0074 ± 0.0006 


0.682 ±0.002 


2.038 ±0.005 


8500 ± 200 


8.446 ±0.079 


24 


01:34:13.8 +30:34:54 


B07141 


2.60 


2.97 ±0.03 


0.026 ± 0.002 


0.004 ± 0.001 


0.348 ± 0.002 


1.042 ±0.005 


8700 ± 800 


8.35 ±0.20 


25 


01:34:13.6 +30:34:49 


B0714g 


2.60 


2.49 ±0.01 


0.0522 ± 0.0008 


0.0053 ± 0.0006 


0.475 ±0.001 


1.409 ±0.003 


8600 ± 300 


8.350 ±0.096 


26 


01:34:13.7 +30:33:45 


S14 


2.82 


3.20 ±0.02 


0.041 ±0.001 


0.0049 ±0.0007 


0.347 ±0.001 


1.046 ±0.003 


9 100 ±400 


8.30 ±0.11 


27 


01:33:17.1 +30:31:11 


C121 


2.94 


2.27 ±0.01 


0.094 ±0.001 


0.0082 ±0.0005 


0.661 ±0.002 


1.955 ±0.005 


8900 ±100 


8.346 ±0.059 


28 


01:34:17.3 +30:33:46 


B0711a 


3.07 


2.54 ±0.02 


0.069 ±0.001 


0.0046 ±0.0009 


0.377 ±0.002 


1.118±0.004 


8800 ±600 


8.28±0.15 


29 


01:34:19.9 +30:33:56 


C192 


3.23 


4. 18 ±0.07 


0.56 ±0.01 


0.053 ±0.003 


1.44 ±0.01 


4.25 ±0.03 


12500 ±300 


8.174±0.059 


30 


01:34:19.3 +30:33:41 


B0712a 


3.24 


3.5 ±0.2 


0.33 ±0.02 


0.038 ±0.007 


1.11 ±0.03 


3.38 ±0.08 


12100 ±800 


8.13±0.13 


31 


01:33:13.9 +30:29:44 


C023 


3.31 


2.44 ±0.02 


0.090 ±0.002 


0.011 ±0.001 


0.635 ±0.002 


1.885 ±0.006 


9700 ±200 


8.224 ±0.067 


32 


01:33:07.6 +30:31:01 


C403 


3.52 


2.26 ±0.02 


0.0720 ±0.0009 


0.01000 ±0.0006 


0.661 ±0.002 


1.970 ±0.005 


9300 ±100 


8.268 ±0.068 


33 


01:33:09.5 +30:29:52 


B0217a 


3.53 


2.35 ±0.03 


0.108 ±0.002 


0.008 ±0.001 


0.783 ±0.004 


2.41 ±0.01 


8400 ±200 


8.48 ±0.10 


34 


01:34:23.2 +30:33:03 


B0715b 


3.64 


4.2 ±0.2 


0.27 ±0.01 


0.031 ±0.008 


0.95 ±0.02 


3.00 ±0.05 


12000 ± 1000 


8.19±0.19 


35 


01:33:14.3 +30:27:11 


C034 


3.68 


2.64 ±0.03 


0.116±0.003 


0.008 ± 0.002 


0.719 ±0.003 


2. 122 ±0.008 


8600 ± 400 


8.44±0.13 


36 


01:33:25.0 +30:25:31 


B0222 


3.70 


3.5±0.1 


0.042 ± 0.003 


0.007 ± 0.002 


0.301 ±0.003 


0.854 ±0.008 


10700 ±900 


8.10±0.18 


37 


01:33:24.7 +30:25:32 


B0222 


3.70 


1.98 ±0.03 


0.060 ±0.003 


0.008 ± 0.002 


0.514 ±0.003 


1.490 ±0.008 


9500 ±700 


8.16±0.15 


38 


01:33:09.8 +30:27:23 


B0221 


3.85 


1.59 ±0.04 


0.163 ±0.005 


0.014 ±0.002 


0.896 ±0.009 


2.71 ±0.03 


9300 ± 400 


8.28±0.10 


39 


01:33:00.5 +30:30:44 


B0218 


4.01 


1.71 ±0.02 


0.144 ±0.003 


0.022 ±0.001 


0.740 ± 0.004 


2.15±0.01 


1 1700 ± 300 


7.926 ±0.064 


40 


01:33:00.3 +30:30:47 


B0218 


4.03 


1.97 ±0.02 


0.1 39 ±0.002 


0.0159 ±0.0008 


0.784 ±0.003 


2.356 ±0.007 


10200 ±200 


8.157 ±0.059 


41 


01:32:57.2 +30:35:57 


B0233c 


4.15 


4.8±0.1 


0.087 ±0.006 


0.013 ±0.004 


0.511 ±0.007 


1.61 ±0.02 


1 1000 ± 1000 


8.23 ±0.19 


42 


01:32:56.9 +30:27:29 


B0243 


4.54 


1.3 ± 0.1 


1.3 ± 0.1 


0.15 ±0.02 


2.8±0.1 


8.0 ±0.3 


14600 ± 700 


8.028 ±0.091 


43 


01:33:09.2 +30:23:30 


B0257 


4.56 


3.31 ±0.04 


0.043 ±0.002 


0.011 ±0.001 


0.610 ±0.003 


1.837 ±0.009 


9800 ±300 


8.290 ±0.090 


44 


01:33:10.3 +30:23:05 


B0255a 


4.61 


3.04 ±0.06 


0.140 ±0.003 


0.022 ±0.001 


0.707 ±0.003 


2.088 ±0.009 


11800 ±200 


8.036 ±0.062 


45 


01:33:07.4 +30:23:13 


B0259a 


4.68 


3.01 ±0.07 


0.027 ±0.004 


0.009 ±0.003 


0.335 ±0.004 


1.045 ±0.009 


11000 ±1000 


8.01 ±0.22 


46 


01:33:06.1 +30:23:27 


B0258b 


4.68 


2.27 ±0.03 


0.464 ±0.007 


0.033 ±0.002 


1.74±0.01 


5.17±0.03 


10000 ±200 


8.418 ±0.072 


47 


01:33:01.6+30:24:31 


C022 


4.69 


1.89±0.01 


0.202 ±0.002 


0.0171 ±0.0006 


1.037 ±0.002 


3.081 ±0.007 


9600 ±100 


8.305 ±0.080 


48 


01:33:45.0 +30:21:38 


B0248 


4.74 


1.300 ±0.007 


0.191 ±0.001 


0.0172 ±0.0005 


1.154 ±0.002 


3.454 ±0.007 


9300 ±80 


8.332 ±0.057 


49 


01:33:00.3 +30:24:33 


C021 


4.75 


2.7±0.1 


0.126 ±0.007 


0.016 ±0.003 


0.86 ±0.01 


2.52 ±0.03 


10000 ±600 


8.26±0.15 


50 


01:32:59.1 +30:24:18 


C018 


4.85 


1.86±0.08 


0.48 ±0.02 


0.052 ± 0.007 


1.86 ±0.03 


5.9±0.1 


11 100 ±500 


8.28±0.10 


51 


01:33:36.6 +30:20:14 


MA19 


4.99 


1.40 ±0.01 


0.227 ±0.002 


0.0209 ± 0.0008 


1.176±0.003 


3.531 ±0.008 


9800 ± 100 


8.273 ±0.062 



TABLE 2 — Continued 



Obj. 




Coordinates" 


Name b 


Rgal 


[onr 


[Ne III] C 


[o ni] c 


[0 III] C 


[o ni] c 


T e ([0 III]) 


12+log([0/H]) 






("2000,(52000) 




(kpc) 


A3727 A 


A3869 A 


A4363 A 


A4959 A 


A5007 A 


(K) 




52 


01 


32:44.4 +30:35:17 


C400 


5.15 


2.78 ±0.03 


0.033 ± 0.002 


0.012 ±0.001 


0.426 ± 0.002 


1.257 ±0.006 


1 1400 ± 400 


7.979 ±0.076 


53 


01 


32:54.6 +30:23:22 


B0261 


5.21 


1.03 ± 0.01 


0.095 ± 0.002 


0.0078 ± 0.0009 


0.810 ±0.004 


2.35 ±0.01 


8300 ± 200 


8.359 ±0.083 


54 


01 


32:43.5 +30:35:17 


C398 


5.22 


4.36 ±0.07 


0.034 ± 0.002 


0.009 ± 0.002 


0.421 ±0.003 


1.263 ±0.007 


10400 ±500 


8.24±0.12 


55 


01 


32:45.9 +30:38:54 


NGC588 


5.23 


0.670 ±0.003 


0.419 ±0.002 


0.0481 ±0.0003 


1.925 ±0.003 


5.663 ±0.008 


10930 ±20 


8.224 ±0.045 


56 


01 


32:45.0 +30:39:02 


NGC588 


5.32 


1.78 ±0.02 


0.215 ±0.003 


0.0191 ±0.0007 


1.17 ±0.01 


3.36 ±0.03 


9600 ± 100 


8.315 ±0.061 


57 


01 


33:30.8 +30:18:38 


B0252 


5.37 


3.00 ±0.03 


0.146 ±0.002 


0.009 ±0.001 


0.851 ±0.003 


2.536 ±0.008 


8600 ± 300 


8.508 ±0.098 


58 


01 


33:29.4 +30:18:06 


B0253 


5.50 


2.72 ±0.09 


0.050 ±0.003 


0.017 ±0.002 


0.688 ±0.007 


2.05 ±0.02 


10800 ±400 


8. 115 ±0.094 


59 


01 


33:21.2+30:17:58 


B0254 


5.56 


2.2 ±0.1 


0.33 ±0.02 


0.03 ±0.01 


1.20 ±0.03 


3.73 ±0.08 


1 1000 ± 1000 


8.16=1=0.19 


60 


01 


32:41.0 +30:24:24 


B0266 


5.83 


3.52±0.05 


0.035 ±0.002 


0.011 ±0.001 


0.522 ±0.003 


1.574 ±0.008 


10300 ±300 


8.216±0.084 


61 


01 


32:39.8 +30:22:28 


C395 


6.12 


2.31 ±0.03 


0.039 ±0.002 


0.013 ±0.001 


0.614 ±0.004 


1.839 ±0.009 


10200 ±300 


8. 129 ±0.085 



Actual position observed based on targets identified in the Local Group Survey images o flMassey et all 12003 ). Closest matching cataloged H II region For observed coordinates. Catalog abbreviations are: MA = |MayalI & AUej 119421 . S = |Searlj 1197 H . B = |Boulesteix et alj |197# , C = |Courtes et alj [19871 . Full designatio 



-ub-lnkiuiv within iin_li\ idual 



H II regions (indicated by the letters at the end of ihe catalog names) are taken from liod.ye el al. 20021. Extinction corrected line fluxes relative to H/3 = 1. The same object was observed on multiple slitmasks. 

TABLE 3 

Balmer Series Line Ratios and Derived Properties 



Obj. 


EWH/3" 

(A) 


c(H/3) 


Stellar Absorption EW 

(A) 




Observed 




Extinction / Absorption 


Corrected 


H7 


H<5 


H7 


H9 


11 ) 


H<5 


H7 


H9 


1 


135.3±0.8 


0.46 ±0.01 


2.0±0.1 


0.390 


0.175 


0.106 


0.032 


0.478 


0.249 


0.171 


0.077 


2 


268 ±2 


0.155 ±0.008 


1.6 ±0.1 


0.436 


0.2312 


0.1200 


0.0552 


0.467 


0.260 


0.142 


0.073 


3 


310±3O 


0.29 ±0.04 


1±1 


0.416 


0.213 


0.239 


0.044 


0.47 


0.26 


0.30 


0.07 


4 


190 ±3 


0.25 ±0.01 


2.9 ±0.2 


0.417 


0.206 


0.099 


0.036 


0.471 


0.256 


0.143 


0.074 


5 


168 ±2 


0.98 ±0.06 


7±1 


0.304 


0.123 


0.079 


0.022 


0.47 


0.26 


0.20 


0.10 


6 


351 =b 3 


0.241 ±0.006 


1.0±0.1 


0.4303 


0.2202 


0.1146 


0.0566 


0.473 


0.255 


0.139 


0.074 


7 


357 ±3 


0.266 ± 0.006 


2.2 ±0.2 


0.4265 


0.2112 


0.1248 


0.0526 


0.476 


0.252 


0.157 


0.075 


8 


640 ±10 


0.02 ±0.05 


8±2 


0.441 


0.229 


0.301 


0.044 


0.47 


0.26 


0.34 


0.07 


8 


440 ±20 


0.00 ±0.02 


2±2 


0.468 


0.255 


0.325 


0.066 


0.471 


0.259 


0.329 


0.072 


Qb 




0.00 ±0.09 


1 1 =b 2 


0.544 


0.285 


0.400 


0.064 


0.52 


0.23 


0.34 


0.004 


10 


193 ±6 


0.43 ±0.05 


2±1 


0.391 


0.193 


0.120 


0.041 


0.47 


0.26 


0.18 


0.08 


11 


110.6±0.5 


0.05 ±0.03 


1.0±0.7 


0.453 


0.202 


0.174 


0.046 


0.47 


0.22 


0.19 


0.06 


12 


750 ± 70 


0.20 ±0.02 


2.2±0.9 


0.428 


0.226 


0.242 


0.055 


0.466 


0.260 


0.283 


0.072 


12 


4800 ± 600 


0.15 ±0.03 


10±2 


0.437 


0.223 


0.241 


0.058 


0.468 


0.26 


0.272 


0.09 


13 


129.5 ±0.3 


0.05 ±0.01 


3.9 ±0.2 


0.426 


0.194 


0.098 


0.026 


0.468 


0.248 


0.155 


0.080 


14 


199 ±2 


1.0±0.1 


3±2 


0.306 


0.144 


0.173 


0.011 


0.47 


0.29 


0.36 


0.09 


15 


203 ±1 


0.220 ±0.008 


0.7 ±0.1 


0.435 


0.2207 


0.1261 


0.0581 


0.474 


0.253 


0.150 


0.075 


16 


116.0±0.4 


0.00 ±0.02 


2.5 ±0.2 


0.456 


0.222 


0.283 


0.042 


0.479 


0.254 


0.316 


0.083 


17 


54.46 ±0.08 


0.00 ±0.01 


3.80 ±0.04 


0.406 


0.162 


0.068 


0.00008 


0.474 


0.256 


0.173 


0.114 


18 


442 ± 1 


0.09 ±0.01 


1.6 ±0.3 


0.452 


0.238 


0.156 


0.063 


0.471 


0.256 


0.172 


0.074 


19 


240 ±10 


0.5±0.1 


1±3 


0.382 


0.142 


0.082 


0.014 


0.47 


0.20 


0.13 


0.03 


20 


248.3 ±0.8 


0.30 ±0.02 


0.9 ±0.6 


0.417 


0.218 


0.228 


0.055 


0.467 


0.260 


0.278 


0.073 


21 


170.6 ±0.3 


0.095 ± 0.004 


2.24 ±0.04 


0.4502 


0.2154 


0.1294 


0.0423 


0.481 


0.249 


0.163 


0.076 


22 


256 ±9 


0.13 ±0.02 


0.4 ±0.3 


0.446 


0.239 


0.131 


0.064 


0.468 


0.259 


0.145 


0.073 


23 


744 ±3 


0.270 ±0.006 


0.0 ±0.3 


0.426 


0.2225 


0.1401 


0.0622 


0.470 


0.257 


0.165 


0.075 


24 


350 ±10 


0.34 ±0.01 


0.0 ±0.2 


0.424 


0.209 


0.113 


0.061 


0.478 


0.249 


0.139 


0.077 


25 


173.4±0.1 


0.288 ±0.007 


0.80 ±0.08 


0.4225 


0.2117 


0.1258 


0.0528 


0.474 


0.254 


0.159 


0.075 


26 


179.7 ±0.2 


0.14±0.01 


1.9=1=0.1 


0.436 


0.219 


0.1083 


0.0438 


0.473 


0.255 


0.140 


0.074 


27 


335 ±4 


0.150 ±0.008 


0.7 ±0.2 


0.446 


0.2328 


0.1450 


0.0631 


0.473 


0.255 


0.162 


0.074 


28 


435 ±8 


0.08 ±0.01 


4.5 ±0.4 


0.446 


0.233 


0.129 


0.051 


0.469 


0.259 


0.154 


0.073 


29 


125.9±0.6 


0.00 ±0.02 


3.1±0.3 


0.450 


0.222 


0.293 


0.046 


0.475 


0.255 


0.33 


0.088 


30 


101 ± 1 


0.08 ±0.05 


1±1 


0.448 


0.205 


0.186 


0.048 


0.47 


0.22 


0.21 


0.06 


31 


191 =b 3 


0.03 ±0.01 


2.5 ±0.2 


0.452 


0.233 


0.145 


0.046 


0.470 


0.257 


0.171 


0.074 
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TAB LE 3 — Continued 



Obj. 


EWH/3 a 


c(H/3) 


Stellar Absorption EW 




Observed 




Extinction / Absorption Corrected 






(A) 




CA1 


H'v 


tio 


rl / 


H9 


rvy 


no 


rl / 


tiq 




32 


250 ±3 


0.211 ±0.009 


2.6±0.1 


0.428 


0.2123 


0.1185 


0.0463 


0.473 


0.254 


0.155 


0.075 




33 


1300 ±200 


0.31 ±0.02 


3.6 ±0.8 


0.420 


0.208 


0.130 


0.052 


0.475 


0.253 


0.167 


0.075 




34 


820 ±70 


0.02 ±0.06 


1±5 


0.462 


0.222 


0.180 


0.073 


0.47 


0.23 


0.19 


0.08 




35 


237 ±1 


0.09 ±0.02 


0.3 ±0.3 


0.452 


0.245 


0.162 


0.067 


0.468 


0.259 


0.173 


0.073 




36 


146.1 ±0.7 


0.15 ±0.04 


2.2±0.8 


0.428 


0.218 


0.117 


0.045 


0.47 


0.26 


0.16 


0.08 




37 


83.2±0.2 


0.00 ±0.02 


3.3 ±0.2 


0.444 


0.203 


0.115 


0.041 


0.479 


0.253 


0.175 


0.106 




38 


600 ±10 


0.90 ±0.03 


0.9 ±4 


0.338 


0.160 


0.110 


0.032 


0.47 


0.26 


0.19 


0.06 




39 


162 ±3 


0.00 ±0.02 


5.7 ±0.2 


0.439 


0.202 


0.128 


0.031 


0.479 


0.254 


0.185 


0.092 




40 


367 ±8 


0.11 ±0.01 


4.4 ±0.3 


0.443 


0.220 


0.145 


0.046 


0.476 


0.252 


0.178 


0.075 




41 


280 ± 30 


0.26 ±0.04 


1±1 


0.421 


0.205 


0.115 


0.025 


0.47 


0.24 


0.14 


0.04 




42 


12.03 ±0.07 


0.5±0.1 


2.5 ±0.2 


0.18 


0.04 


-0.27 


-0.33 


0.47 


0.40 


0.09 


0.0003 




43 


218±7 


0.10±0.02 


2.0 ±0.2 


0.448 


0.228 


0.140 


0.052 


0.474 


0.254 


0.166 


0.075 




44 


51.6±0.2 


0.00 ±0.03 


2.56 ±0.06 


0.437 


0.182 


0.104 


0.006 


0.485 


0.251 


0.183 


0.097 




45 


117±4 


0.00 ±0.03 


4.3 ±0.3 


0.445 


0.191 


0.117 


0.006 


0.49 


0.250 


0.179 


0.083 




46 


137 ±3 


0.00 ±0.02 


7.0 ±0.4 


0.425 


0.193 


0.244 


0.037 


0.469 


0.259 


0.301 


0.109 




47 


250 ±2 


0.054 ±0.008 


1.6±0.1 


0.456 


0.2369 


0.1896 


0.0558 


0.473 


0.255 


0.209 


0.074 




48 


216.2 ±0.2 


0.155 ±0.007 


2.2±0.1 


0.4395 


0.2167 


0.1646 


0.0490 


0.476 


0.252 


0.200 


0.075 




49 


119±4 


0.29 ±0.06 


2.4 ±0.8 


0.401 


0.195 


0.136 


0.032 


0.47 


0.26 


0.20 


0.08 




50 


27.6 ±0.4 


0.25 ±0.05 


1.5 ±0.2 


0.373 


0.124 


0.000 


-0.045 


0.47 


0.23 


0.10 


0.06 




51 


214.9 ±0.3 


0.273 ±0.009 


3.0±0.1 


0.420 


0.1923 


0.1659 


0.0384 


0.481 


0.246 


0.222 


0.077 




52 


138 ±2 


0.07 ±0.02 


2.7 ±0.2 


0.445 


0.213 


0.094 


0.040 


0.478 


0.249 


0.131 


0.077 




53 


121 ±1 


0.08 ±0.02 


6.3 ±0.2 


0.417 


0.167 


0.074 


0.028 


0.478 


0.247 


0.161 


0.085 


54 


250 ±10 


0.18±0.02 


3.2 ±0.4 


0.432 


0.211 


0.101 


0.044 


0.475 


0.252 


0.136 


0.076 


in 


55 


455 ±3 


0.082 ±0.005 


1.83 ±0.07 


0.4539 


0.2316 


0.2481 


0.0558 


0.476 


0.253 


0.276 


0.074 


o" 


56 


3 100 ±500 


0.418 ±0.009 


10 ±10 


0.421 


0.184 


0.156 


0.0483 


0.50 


0.24 


0.21 


0.07 




57 


243.3 ±0.7 


0.18±0.01 


0.4 ±0.2 


0.439 


0.231 


0.173 


0.062 


0.470 


0.257 


0.196 


0.074 


:/> 

■3* 


58 


117.3±0.3 


0.09 ±0.04 


5.3 ±0.7 


0.415 


0.194 


0.120 


0.034 


0.47 


0.26 


0.19 


0.11 


59 


20.76 ±0.04 


0.22 ±0.07 


1.1 ±0.2 


0.37 


0.088 


0.00 


-0.08 


0.47 


0.19 


0.11 


0.04 




60 


164 ±4 


0.23 ±0.02 


2.4 ±0.2 


0.422 


0.202 


0.105 


0.039 


0.475 


0.252 


0.148 


0.075 


61 


95 ±1 


0.00 ±0.01 


4.9 ±0.1 


0.420 


0.186 


0.107 


0.024 


0.473 


0.257 


0.188 


0.112 


§ 



Note. — Line fluxes are given relative to H/3=l and rounded to the precision of one significant figure in the error. 
a The equivalent width of H(3 in emission.* 3 Because there is no discernible continuum emission associated with the nebular line emission from this H II region, the equivalent width of the H/3 emission is undefined. 



